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Abstract: We present a new technique for computing Hilbert series of A/" = 1 super- 
symmetric QCD in four dimensions with unitary and special unitary gauge groups. We 
show that the Hilbert series of this theory can be written in terms of determinants of 
Toeplitz matrices. Applying related theorems from random matrix theory, we compute 
a number of exact Hilbert series as well as asymptotic formulae for large numbers of 
colours and flavours - many of which have not been derived before. 
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1 Introduction and summary 

Given a supersymmetric gauge theory, the space of solutions of the vacuum equations, 
known as the moduli space of vacua or simply the moduli space, is one of the first 
properties one can study. The phase structure and the possible excitations at a given 
vacuum configuration can be investigated. In supersymmetric gauge theories, holo- 
morphic gauge invariant operators play a central role in giving us information about 
the structure of the moduli space. It is therefore worthwhile to study and count these 
holomorphic functions. This can be done in a similar fashion to statistical mechanics, 
i.e. by constructing a partition function. Mathematically, this partition function is 
known as the Hilbert series (see e.g. [1, 2] for reviews of applications of the Hilbert 
series to gauge theories and string theory). The Hilbert series not only contains infor- 
mation about the spectrum of the operators in the theory, it also carries geometrical 
properties of the moduli space (see e.g. [3, 4, 9]). It is also an indicator of whether the 
moduli space is Calabi-Yau (see e.g. [5, 9, 10]). Moreover, Hilbert series can be used 
as a primary tool to test various dualities in gauge theories and in string theory (see 
e.g. [6, 8, 17, 24]). 

The Hilbert series was first applied to study the moduli space of M = 1 Super- 
symmetric Quantum Chromodynamics (SQCD) in four dimensions by Pouliot in 1998 
[7]. In that paper, a method for computing the Hilbert series using the Molien-Weyl 
formula was briefly discussed and the Hilbert series of SU(2) SQCD up to 3 flavours 
were explicitly stated. In 2005, Romelsberger [8] computed the Hilbert series for SU{2) 
with 3 flavours and used it to verify the Seiberg duality. Later, in [9, 10], the Hilbert 
series were explicitly computed for SQCD with classical gauge groups for various num- 
bers of colours and flavours and from which several geometrical aspects of the moduli 
space were extracted. In addition to SQCD, there have been several applications of 
the Hilbert series to study a wide range of supersymmetric gauge theories and string 
theory, e.g. quiver gauge theories on D3-branes [11, 12], Af = 4 gauge theory in four di- 
mensions [13], SQCD with one adjoint chiral superfield [14], M2-brane theories [15-19], 
perturbative string spectrum [20], moduli spaces of instantons [21, 22], and a number 
of Af = 2 gauge theories in four dimensions [22-24] . 
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In this paper, we focus on the Hilbert series of SQCD with the gauge groups U(N C ) 
and SU(N C ) and Nf flavours of fundamental matters. Similarly to [9], the Molien- 
Weyl formula is used as a starting point of the computations. However, instead of 
evaluating the Hilbert series case by case for each (Nf,N c ) as in [9], we take another 
approach of computations. First, we show that the Molien-Weyl formula for U(N C ) and 
SU(N C ) SQCD can be rewritten in terms of determinants of Toeplitz matrices. Then, 
we evaluate these determinants both exactly and asymptotically for a large class of 
values of Nf and N c using several theorems from random matrix theory. 1 In particular, 
such theorems are stated in §2.2.2 for U(N C ) SQCD and in §3.2 for SU(N C ) SQCD (see 
e.g. [40] for a mathematical review). 

Since a different method from [9] is used to compute Hilbert series in this paper, the 
way our results presented here is generally distinct from that in [9]. In this paper, the 
Hilbert series are computed for a class of values of Nf and N c rather than a specific value 
of (Nf,N c ) - for example, the exact Hilbert series for SU(N C ) SQCD with Nf = N c 
and Nf = N c + 1 are presented respectively in §3.2.3 and §3.2.4 and such results, of 
course, apply for any iV c > 1. As a consequence, many general results conjectured 
in [9] can be proven (or at least can be checked in a non-trivial way) using Toeplitz 
determinants. Such a technique also allows us to compute explicity asymptotic formulae 
for large Nf and N c . A number of exact and asymptotic results presented in this paper 
are too difficult or impossible to be derived using the method in [9], partly due to a 
large number of contour integrals in the Molien-Weyl formula and partly due to the 
complications of the results themselves. 

Let us summarise the main results and key features of this paper below. 

Outline and key Points: 

• In §2.2.1, we give the definitions of a Toeplitz matrix and a Toeplitz determinant. 
In §2.2.2, we state important facts, namely the CGBO formula and the strong 
Szego limit theorem, which are applied to compute Toeplitz determinants for 
U(N C ) SQCD with N f flavours. 

• In §2.2.3, we show that the Hilbert series for U(N C ) SQCD with Nf flavours can 
be written in terms of a Toeplitz determinant of a symbol with a zero winding 
number. 

• In §2.4, we use the CGBO formula to compute exact Hilbert series for U(N C ) 
SQCD with N f < N c , N f = N c + 1, N f = N c + 2 and N f = N c + 3. We also 

1 We mention en passant that there have been several recent works on applications of matrix models 
to supersymmetric gauge theories and string theory - see, for example, [50-56]. 
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conjectured the general expression (2.79) of the Hilbert series for any N c and any 
Nf written in terms of a sum over representations of SU(Nf) x SU(Nf). The 
information about the moduli space and the chiral ring of U(N C ) SQCD can be 
extracted from these Hilbert series - this is summarised in §2.1. 

• In §2.6.1, we study asymptotics for Nf,N c » 1 with a fixed finite difference 
Nf — N c > 0. The asymptotic formula for this limit is given by (2.94). In 
§2.6.2, we study asymptotics for Nf, N c » 1 with a fixed ratio Nf/N c > 1. The 
asymptotic formula for this limit is given by (2.104). 

• In §3.1, we show that the Hilbert series for SU (N c ) SQCD with Nf flavours can be 
written in terms of a sum of 3 parts: (i) the Hilbert series of U(N C ) SQCD with Nf 
flavours, (ii) Toeplitz determinants of symbols with positive winding numbers, and 
(iii) Toeplitz determinants of symbols with negative winding numbers. Observe 
that we can use previous results from U(N C ) SQCD in part (i) of the SU(N C ) 
SQCD Hilbert series computations. 

• In §3.2, we state important facts, namely the Bottcher-Widom theorem and the 
Fisher-Hartwig theorem, which are applied to compute Toeplitz determinants for 
SU(N C ) SQCD with N f flavours. 

• In §§3.2.2, 3.2.3 and 3.2.4, we use the Bottcher-Widom formula to compute exact 
Hilbert series for SU(N C ) SQCD with N f <N c - 1, N f = N c , N f = N c + 1. The 
first two subsections contain the proofs of (3.6), (5.3), (3.15) and (5.2) of [9]. 
These exact results also provide a non-trivial test for the general formula (3.58). 

• In §3.4, we apply the Fisher-Hartwig Theorem to examine asymptotics of the 
Hilbert series of SU(N C ) with Nf flavours when Nf,N c » 1. In §3.4.1, we 
study asymptotics for Nf,N c » 1 with a fixed finite difference Nf — N c > 0. 
The asymptotic formula for this limit is given by (3.67). In §3.4.2, we study 
asymptotics for Nf,N c » 1 with a fixed ratio Nf/N c > 1. The asymptotic 
formula for this limit is given by (3.74). 

Note added. In the second version of this paper, the authors become aware of the use 
of Toeplitz determinants and related theorems in the literature on decaying D-branes; 
see e.g. [57-61]. 

Notation for representations and characters. 

We denote an irreducible representation of a group SU(r + 1) by a Dynkin label 
[d, a 2 , ■ ■ ■ , a r ]. For example, we denote by [1,0,..., 0] the fundamental representation 
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and by [0, ...,0,1] the anti-fundamental representation. We also use the subscripts 
k; L and k; R to indicate respectively the k-th postitions from the left and the right, 
e.g. Ik-L in [0, . . . , 0, lk-L, 0, . . . , 0] denotes the 1 in the fc-th position from the left. 

For representations of the product group Gi x G2, we use the notation [...;...] 
where the tuple to the left of the ; is the representation of Gq, and the tuple to the 
right of the ; is the representation of G 2 . 

Since a representation is determined by its character, we denote by [ai, a 2 , ■ ■ ■ , a r ] x a 
character of the representation [ai, ei2, . . . , a r ] written in terms of the variable x±, . . . , x r . 
Here the subscript x denotes collectively the set of variables aq, . . . , x r . For example, the 
characters of the fundamental representation and the anti-fundamental representations 
of SU(Nf) are given by (2.10) and (2.11). 

2 U(N C ) SQCD with N f flavours 

Consider Af — 1 supersymmetric U (N c ) gauge theory in four dimensions with Nf quark 
Q\ and Nf anti-quarks Q", where a = 1, . . . , N c and i = 1, . . . , Nf. Let us take the 
superpotential of this theory to be zero: W = 0. The information about the gauge 
and global symmetries as well as how the matters transform under such symmetries is 
collected in Table 1 (see, e.g., [33]). 





Gauge symmetry 
U(N C ) = SU(N C ) x 17(1) 


Global symmetry 

SUiNfh SU(N f ) 2 U(1) B U(1)r U(1)q U(1)q 


Qa 

02 


[0,... ,0,15-1] 
[1,0,...,0;+1] 


[1,0,. ..,0] [0,...,0] 1 ^ 1 
[0,...,0] [0,...,0,1] -1 -1 




Table 1. The gauge and global symmetries of SQCD and the quantum numbers of the chiral 
supermultiplets. The quarks are Q\ while the antiquarks are Q^. We also draw it as a quiver 
theory. The circular node represents the U(N C ) gauge symmetry while the two square nodes 
represent global U(Nf)\ and U(Nf) 2 symmetries. Each square node gives rise to a baryonic 
J7(l) global symmetry, one of which is redundant. We thus have U(1)q q that combine into 
the non-anomalous U(1)b (sum) and anomalous U{1)a (difference). 

ftf — 1 SQCD with U (N c ) gauge group posesses several interesting phenomena, 
e.g. the Seiberg duality [33], non-Abelian flux tubes (strings) and confinements [33, 34]. 
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Another motivation for studying this theory is that, as we shall see in §3.1, Hilbert series 
for U(N C ) SQCD with Nf flavours are needed for our computations of Hilbert series of 
SU(N C ) SQCD with Nf flavours. 

2.1 The moduli space of vacua 

In this paper, we focus on the classical moduli space of U(N C ) SQCD with Nf flavours. 

The case of Nf < N c . At the generic point of the moduli space, the U(N C ) gauge 
symmetry is partially broken to U(N C — Nf). Thus, there are 

N 2 C - (iV c — Nf) 2 = 2N c N f - N 2 f 

broken generators. The total number of degrees of freedom of the system is, of course, 
unaffected by this spontaneous symmetry breaking and the massive gauge bosons each 
eat one degree of freedom from the chiral matter via the Higgs effect. Therefore, of the 
original 2N c Nf chiral multiplets (quarks and antiquarks), only 

2N c N f - (2N c N f — Nj) — N] 

gauge invariant degrees of freedom are left massless. Hence, the dimension of the moduli 
space is 

M Nf<Nc = N 2 . (2.1) 

We describe these Nf massless degrees of freedom by a collection of gauge invariant 
objects, called mesons: 

M) = Q\Q) , (2.2) 
where the mesons transform in the bi-fundamental representation 

[1,0,...,0;0,...,0,1] 

of SU(Nf)i x SU(Nf) 2 . The moduli space is said to be freely generated by the 
mesons. 

The case of Nf > N c . At the generic point of the moduli space, the U(N C ) gauge 
symmetry is completely broken. Thus, there are N 2 broken generators. Therefore, of 
the original 2N c Nf chiral multiplets (quarks and antiquarks), only 

2N c N f - N 2 
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gauge invariant degrees of freedom are left massless. Hence, the dimension of the moduli 
space is 

dim M Nf > Nc = 2N c N f - N 2 C . (2.3) 

The moduli space is still generated by the mesons: 

Mj = Q\Q°- . (2.4) 

However, for Nf > N c + 1, there are non-trivial relations between these mesons. The 
basic relations (which generate the entire ideal of relations) can be written as 

. . . = . (2.5) 

They transform in the representation 

[0, . . . , 0, l Nc+1 , L , 0, . . . , 0; 0, . . . , 0, l^+i;*, 0, . . . , 0] (2.6) 

of the SU(Nf)i x SU(Nf)2 flavour symmetry. Note that the dimension of this repre- 
sentation is the number of relations between the mesons. 

The case of Nf = N c + 1. In this case, the representation (2.6) reduces to the trivial 
representation, i.e. we have precisely one relation. Equation (2.5) becomes 

detM = 0. (2.7) 

From (2.3), the dimension of moduli space is 

Nf — 1 . (2.8) 

Since this is the number of generators (which is Nj) minus the number of relations 
(which is 1), the moduli space for the case Nf = N c + 1 is a complete intersection. 

2.2 The computations of Hilbert series 

The Hilbert series of U(N C ) SQCD with Nf flavours can be computed in two steps as 
follows. 

Step 1. First we consider the space of symmetric functions of quarks Q and anti- 
quarks Q. The Hilbert series of this space can be constructed using the plethystic 
exponential, which is a generator for symmetrisation [3, 4]. To remind the reader, we 
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define the plethystic exponential of a multi-variable function g{t\, ...,t n ) that vanishes 
at the origin, g(Q, 0) = 0, to be 



PEfofa, . . . , t n )} := exp (j2 • • • > tn) J 



(2.9) 



Let t be the U (1)q global charge fugacity and t be the U(1)q global charge fugacity. 
Note that < t, t < 1. Let z%, . . . , zn c be the fugacities of the U(N C ) gauge group. The 
character of the fundamental representation [1, 0, . . . , 0]+i of U (N c ) can be written as 
J2a=i z a> whereas the character of the antifundamental representation [0, . . . , 0, l]_i of 
U (N c ) can be written as Ylia=i z a l 

Let Xi, . . . , xn, be the SU(Nf)i fugacities and let yx, ■ ■ ■ , Un s be the SU(Nf) 2 fugac- 
ities. Explicitly, we take the character of the fundamental representation of SU(Nf)i 
to be 



N,-2 



[l,0,...,0] 1 . = x 1 + J2 



1 



, =1 %k x Nf—l 

and take the character of the anti-fundamental representation of SU(Nf) 2 to be 

N f -2 



(2.10) 



[o,o,...,i] y = i+ £ + y Nf -i 
yi t^i Vk+l 



(2.11) 



Then, the Hilbert series of the space of symmetric functions of quarks Q and 
antiquarks Q can be written as 



PE 



N c 



N c 



0=1 



a=l 



(2.12) 



where t[l,0, . . . ,0} x Y^ a =i z a 1 comes from the quarks Q and t[0, 0, . . . , l] y 5^o=i z <* comes 
from the antiquarks Q. Using the definition of the plethystic exponential and the 
identity — log(l — x) — J2T=i xk /k' we can wr ite down the expression (2.12) explicitly 
as a rational function: 



N c 

n 



i 



t\ (i - tz?x x ) nfir 2 (i - tz-^) (i - te-^^-i) 

N c 1 

TT L 

:i (i - t Za A nfir 2 (i - ^) (i - 



X 



(2.13) 
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One may set xi, . . . , Xn c , yi, . . . , y^ c to unity (this process is called an unrefinement) 
and obtain 



N c 1 

U {1 _ i)Nf{1 _ ha)Nf ■ ( 2 - 14 ) 

Step 2. Since the moduli space is parametrised by gauge invariant quantities, we 
need to project representations associated with symmetric functions in Q and Q dis- 
cussed in Step 1 onto the trivial subrepresentation, which consists of the quantities 
invariant under the action of the gauge group. Using knowledge from representation 
theory (known as the Molien-Weyl formula - see e.g. [35, 36]), this can be done by 
integrating over the whole gauge group. In this section, we are interested in the U(N C ) 
gauge group, whose Haar measure is given by (see e.g. [37, 38]) 

dzx f dz Nc , A 2 



N c \[2m) N - J\ Zl \ =1 z x J\z Nc \=i z Nc 
where 

\A Nc (z)\ 2 = J] \z a -z h \ 2 . (2.16) 

l<a<b<Af c 

The Hilbert series for U(N C ) SQCD with Nf flavours is given by 



9N f ,u(N c ){t,t,x,y) = J dfx u{Nc) PE 



N c N L 

Z, 



t[l,0,...,0} x Y,^a 1 +t[0,0,...,l] y J2 

a=l a=l 



(2.17) 



where here and henceforth we write x and y as collective notation for x%, . . . , x^ c and 
yi, ■ ■ ■ i Vn c - Setting x\, . . . , xn c and y\, . . . , yN c to unity, we obtain the unrefined Hilbert 
series of U(N C ) SQCD with Nf flavours: 

f Nc 1 

9N flU (N c )(t,t) = / d mNc) J] ± N ~ Nf ■ (2-18) 

a=l v z a ' \ a ' 

2.2.1 Toeplitz matrices and Toeplitz determinants 

In this paper, we are concerned with computing the integrals (2.17) and (2.18) by means 
of Toeplitz matrices and their determinants. In this section, let us briefly summarise 
important facts which will be useful in subsequent computations. We follow [40] closely 
and neglect subtleties regarding convergences. We also refer the reader to [39] for a 
comprehensive review on Toeplitz operators. 
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Let T be a complex unit circle and let : T — > C be a continuous function such 
that <p(z) for all z G T and has a zero winding number around zero, i.e. 

±m. (2 . 19) 

2mJ ]z]=1 z 4>{z) 
The Fourier coefficients <pk (with k G Z) of / is defined by 



d6f(e l »)e- lkti = — f —z-*f{z). (2.20) 
2tt ./o 2tu J n=1 2 



The n x n Toeplitz matrix is defined by 

T n (0) := {^-k)l k=l ■ (2.21) 

The function is sometimes referred to as the symbol of the Toeplitz matrix T n . We 
are also interested in computing the determinant of the Toeplitz matrix, 

D n (4>) :=detT n (0) . (2.22) 

This is called the Toeplitz determinant of the symbol 0. We also define the infinite 
Toeplitz matrix T(0) and the infinite Hankel matrix H(<p) as 

W) := , H{<f>) := (^■ +Jfc _i)~ fcsl • (2.23) 

Subsequently, we find that the following definition is useful 

H(0) := (0_,_ fc+1 )^ =1 . (2.24) 

2.2.2 The Geronimo-Case-Borodin-Okounkov (GCBO) formula and the 
strong Szego limit theorem 

There are theorems which are very useful for computing Toeplitz determinants (see, 
e.g., [40, 41]). Before stating these theorems, let us give more definitions which will be 
extensively referred to later. 

Under general conditions can be written as 

= 0+0- , (2.25) 

where + (resp. 0_) is a nonzero analytic function in the interior (resp. exterior) of T. 
This is factorisation is known as the Wiener-Hopf factorisation. Let us define the 
functions 

U := 5^ , V := . (2.26) 
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Define the Fourier coefficients of U and V as 
We also define 

G{4>) := exp(log 4>) = exp j ^ dz z~ l log ^ . (2.28) 

Define the projection matrices P k and Q k to be infinite matrices such that 

P k = diag( l,l, . ..,l , 0, 0, . . .) , Q k = diag(0,0 ,0, 1, 1, . . .) . (2.29) 

k ones fc zeros 

Now let us state the theorems. For the proofs, we refer the reader to [40]. 

Theorem 2.1. (The Geronimo-Case-Borodin-Okounkov (GCBO) formula. 
[43, 44]) The operator 1 — H(b)H(c) is invertible and the n x n Toeplitz determi- 
nant of the symbol <fi is given by 

A>(0) = G(0) n £(0)det(l - Q n H(b)H(c)Q n ) 

= G((f>) n E(<p)desb(l-K n ) , (2.30) 

where the (i,j)-entry of the infinite matrix K n is 

K„(iJ) = {P' a,+ ^- fc,/ ^-" (2.31) 
^ otherwise , 

and 

E(0) = exp^A;(log0) fc (log0)_ fc j . (2.32) 

Theorem 2.2. (The Szego strong limit theorem. [45]) In the limit of large n, 
we have 

D n {<j>) ~ G{<t>) n E{<j>) , (2.33) 
where G(<f>) is defined as in (2.28) and E(<p) is defined as in (2.32) . 

In the following subsections, we give explicit examples on how to apply these the- 
orems to compute Hilbert series of U(N C ) SQCD with Nf flavours. 
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2.2.3 The Hilbert series and the Toeplitz determinant 

Now let us go back to our problem of computing (2.17) and (2.18). Let 

<t>{t,t,x,y,z) = PE[t[l,0,...,0} x z~ 1 +t[0,0,...,l} y z] . (2.34) 

Note that <f>(t, t, x, y, z) is non-zero for all z G T and it has a zero winding number 
around z = 0. The Hilbert series can be written as 

* N c 

9N f ,u(N c )(t,t,x,y) = / dti U ( Ne )(z 1 ,...,z»r c )Y[<l>{t,t,x,y,z a ) . (2.35) 

J a=l 

It is well-known that the U(N C ) Haar measure can be written in terms of a determinant. 
For the sake of completeness and ease of reading, we include a proof in Appendix A. 
In particular, from (A. 6), we find that 

d mNc) = det ( -L / — z a - b ) , (2.36) 



y 2m 

and, from (A. 9), we have 



\ z \ =1 ' Ka,b<N, 



g Nf ,u(N c )(t,t,x,y) = det ( ^— I —z a b <p(t,t,x,y, z) 

\ 2m J\z\=l Z Sl<a,b<N c 

= D N M • (2.37) 

This is simply the N c x N c Toeplitz determinant with the symbol <j)(t,t,x,y, z). We 
therefore take n in Theorem 2.1 and Theorem 2.2 to be the number of colours iV c . 

Results for any Nf and any N c . Before we proceed, let us state some useful results 
for future references. The following statements are true for any Nf and any iV c , 

G((f>) = exp ( -L I — log<f)(t,t,x,y,z)) = 1 , (2.38) 

E{(j>) = PE [[1, 0, . . . , 0; 0, . . . , 0, l] x , y tt\ , (2.39) 

where the first equality can be easily checked and the second equality will be proven in 
the next subsection. 

Let us now look at various cases of Nf and N c . 
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2.3 The strong Szego limit theorem for the iV c >> Nf limit 

In the case, it follows from the strong Szego limit theorem that 

oo 

k>E9N »N f (t,t,x,y) ~ ^ A;(log0) fc (log0)_ fc , (2.40) 

fe=i 

where, for n G Z, (log0) n is defined to be the n-th Fourier coefficient of log</>: 

(log <P) n :=^-.l dzz-"- 1 log <f) . (2.41) 

2m J\z\=l 

This can be evaluated exactly as 

f i[0,...,0,l]^ if n>0 
(log0) n =^O ifn = (2.42) 

[ ^[l,0,...,0] xW tW ifn<0. 

where [0, . . . , 0, l] y n = ± + ^ + ... + and similarly for [1,0,..., 0] x \„\ . 

Hence, it follows from (2.32) that 

E{(j>) = exp ^ ^ [1, 0, . . . , 0; 0, . . . , 0, 1]^* {ttf^j 

= PE[[l,0,...,0;0,...,0,l] aw tt] . (2.43) 
It follows immediately from Theorem 2.2 and (2.38) that 

9N c » Nf (tJ, x, y) ~ E{(f>) = PE [[1, 0, . . . , 0; 0, . . . , 0, l] x;y tt\ . (2.44) 
This formula can also be written in terms of sums of irreducible representations as 

/ ~ 1 
9N c »N t {t,t,x,y) ~ - _ -~— x 

22 K^2,- • • ; n Nf _i,. . . ,n 2 ,n 1 ](tt) T '^ 1 3U] . (2.45) 



ni,...,n ]V/ _i=0 



The unrefined Hilbert series. Setting x\,..., xn c and yi, ■ ■ ■ yN c to unity in (2.43), 
we replace the representations by their dimensions. Note that 

dim[l, 0, . . . , 0] = dim[0, ... ,0,1] = Nf . (2.46) 
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Hence, we have the unrefined Hilbert series: 

9N c »N f (t,t, 1, 1) ~ PE[N 2 f tt] = ?— -s . (2.47) 

(1 - tt) N f 

In the next subsection, we use the GCBO formula to show that (2.43) and (2.47) 
are actually exact for Nf < N c . In other words, corrections to the asymptotic 
formula (2.43) and (2.47) are indeed zero for Nf < N c . 

2.4 The GCBO formula and exact results 

We perform the Wiener-Hopf factorisation of <fi as 

(j)(t,t,x,y,z) = <f> + (t,y,z)</>-(t,x,z) , (2.48) 

where 

Mt, y, z) = pe \t[0, 0, . . . , l\ v z] , 0_ (t, x, z) = PE [t[l, 0, . . . , O]^" 1 ] . 

(2.49) 

Using (2.27), we find that U n and V n can be evaluated exactly. For < n < Nf, 

N f -n 

U n = (-t) n j2 [m, 0, . . . , ; 0, . . . , 0, l (n+m)]R , 0, . . . , Q] xiV {-tt) m , 

m=0 

Nj-n 

V. n = {-t) n [0, • • • , 0, l ( n +m) ;L, 0, . . . , ; 0, . . . , 0, m} x , y {-tt) m , (2.50) 

where the subscripts (n + m); L and (n + m); R indicate respectively the (n + m)-th 
position from the left and the right. For n > Nf + 1, we have 

U n = V. n = . (2.51) 

Unrefined results. Setting xi, . . . ,x^ c and yi, ■ ■ ■ yN c to unity in (2.50), we replace 
the representations by their dimensions. Note that 

Nf + m—1 
m 



dim[m, 0, . . . , 0] = dim[0, . . . , 0, m] - 
dim[0, ... ,0, l (n+m ) ;L ,0, ... ,0] = dim[0, . . . , 0, l( n+m);jR ,0, . . . ,0] = yj^^ ■ 



(2.52) 
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For < n < Nf, they can be written in terms of hypergeometric functions: 

N f -n 



z — 4 V rn I \n + ml 



-tty 



t) n yj\Fi{n-N f ,Nf-,n + l-,tt) , (2.53) 

Nf-n 

«)" E ( J 

to— n \ 



= (-t) n ( N ^2F 1 (n- N f ,N f ;n + l;tF) . (2.54) 
For n > Nf + 1, we have 

U n = V. n = . (2.55) 

The exact Hilbert series. The GCBO formula states that the exact formula for 
9N f ,u(N c ), for any N c and Nf, is given by 

9N f ,u{N c )(t,t, x, y) = det(l - K Nc ) PE [[1, 0, . . . , 0; 0, . . . , 0, l] X]y tt\ . (2.56) 

Hence, the Hilbert series for any iV c and Nf is simply the Hilbert series for Nf < N c 
multiplied by a correction factor det(l—K Nc ). In the following subsections, we compute 
this correction exactly and asymptotically in various cases. 

2.4.1 The case of N f < N c 

Let us apply the GCBO formula to the case of Nf < N c . From (2.31) and (2.51), we 
have 

K Nc {i,j) = foralU,j, (2.57) 

It is easy to see that 

det(l - K Nc ) = 1 . (2.58) 
Thus, from (2.56), the refined Hilbert series is given by 
9N f < Nc {t, % x, y) = PE [[1, 0, . . . , 0; 0, . . . , 0, l] x . y tt\ 

= - 72K^ z2 [ni,n 2j ... J n Nf _ l ;n Nf _ u ... J n 2j n 1 }(tt) T '^ ^ . (2.59) 

v I ni,...,njvx-i=0 
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The unrefined Hilbert series is 

9N f <N c (t, t, 1, 1, 



'I -tip 



(2.60) 



Indeed, for Nf < N c , the moduli space is of Nj dimensional and is freely generated by 



the mesons. 



2.4.2 The case of N f = N c + 1 

Let us focus on the case of Nf = N c + 1. Using (2.31), (2.50) and (2.51), we find that 



Un c +iV-n -i 




if i = N c and j = N c 
otherwise . 



Thus, we have 



t N n N f m = j = N c 

otherwise , 



and so 



det(l - K Nc ) = 1 - (tt) 



Nr 



(2.61) 



(2.62) 



Thus, from (2.56), the refined Hilbert series is given by 

9N f ,u(N f -i)(t, t, x, y) = (l- {tiff) PE [[1, 0, . . . , 0; 0, . . . , 0, l] x . y it\ 

= 22 [ni,n 2} ... } n Nf - 1 ;n Nf - 1} ... } n 2} n 1 }(tt)^ : ^ 1 ]nj . (2.63) 

ni,... ,71^-1=0 

Setting x's and y's to unity, we see that the unrefined Hilbert series is 

_ 1 _ (tt) N f 

g Nf ,u(N f -i)(t,t,l,l) = _ v J 2 . (2.64) 

2.4.3 The case of N f = N c + 2 

Let us focus on the case of Nf = N c + 2. Using Eq. (2.31) and (2.51), we find that 

Un c +\V^ Nc -i + U Nc+2 V- Nc -2 if i = j = N c 
U Nc+1 V- Nc - 2 if i = N c , j = N c + 1 

K Nc (hj) = { U Nc+2 V- Nc -! if % = iV c + 1, j = N c 

U Nb+2 V- Nb -2 if i = N c + 1, j = N c + 1 

otherwise . 
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N f +1 



N, 



Using (2.50), we find that 

K Nc (N c , N c ) = K Nc (N f -2,N f - 2) 

= [i,o,...,ouo,...,o,rj y (tt) 

- ([1, 0, . . . , 0, l] x + [1, 0, . . . , 0, 1]„ + 1) (tQ 
+ [0,...,0,l] I [l,0,...,0] J/ (tt)^- 1 , 
K Nc (N c ,N c + l) = K Nc (N f ,N f -l) 

= [i,o,..., o] x i N ft N f +1 -[i,o,..., o^-H^ 

K Nc (N c + 1,N C )= K Nc (N f - 1, N f ) 

= [0, 0, . . . , l]yt N ^ + l - [0, 0, ... , llt^-^f 

K Nc (N c + 1,N C + 1) = K Nc (N f — l,Nf — 1) 
= (tt) Nf , 

and K^ c {i,j) = 0, otherwise. Therefore, we obtain 



(2.65) 



det(l - K Nc ) = 1 - [0, . . . , 0, 1; 1, 0, . . . , 0] x . y (tt) 



Nf-l 



+ ([l,0,...,0,l;0,...,0U+[0,...,0;l,0,...,0,l] :c . 2/ ) (ttff 
-[l,0,...,0;0,...,0,l]^(tt)^ +1 + (tt) 2 ^ , (2.66) 

and the exact Hilbert series for Nf = N c + 2 is 

9N f ,u(N f -2)( l t,t,x,y) = PE [[1, 0, . . . , 0; 0, . . . , 0, l\ x . y tt\ x 

1 - [0, . . . , 0, 1; 1, 0, . . . , Oj^Ctt)^- 1 

+ ([l,0,...,0,l;0,...,0] a;;?/ + [0,...,0;l,0,...,0,l]^) {tt) N f 

-[l,0,...,0;0,... J 0,l] se . 1 ,(tt) JV ' +1 + (tt) 2JV '] . (2.67) 

This can be rewritten in terms of a sum of representations of SU(Nf) x SU(Nf) as 



9Nf,u(Nf-2){t,t,x,y) = ^2 ( tt ) Eo 



X 



ni,n2,...,n Nf -2=0 



[m, n 2 , . ■ ■ , n Nf - 2 , 0; 0, n N} - 2 , ■ ■ ■ , n 2 , ni] x;y . 
The unrefined Hilbert series. Setting x's and y's to unity, we have 
' (tt) N f + Nj{tt) N f- l {\ - ttf iii=j = N f -2 



(2.68) 



K Nc {hj) = < 



-N f t- l {tt) N f{l-tt) 
-N f t- l {tt) N s(l-tt) 



(tt) 




N f 



if % = N f - 2,j = N f -1 

if % = N f - l,j = N f - 2 (2.69) 

if % = j = N f - 1 

otherwise . 
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Hence, we arrive at 

det(l - K Nc 



(l - (tt) N, y - N 2 f (tt)^- 1 (1 -tt) 2 . (2.70) 



The unrefined Hilbert series for the Nf = N c + 2 is therefore 



9N f ,U(N f ~2)(t, t) 



(1 - tt) N f 

Examples. We list a few examples of unrefined Hilbert series for small Nf. 

1 + AU + t 2 t 2 



(2.71) 



93, U(l)(t,t) 

94, U{2)(t,t) 

95, U(3)(t,t) 

96, U(4)(t,t) 



(i-ttf • 

1 + Att + 10t 2 t 2 + At 3 t 3 + t 4 t 4 



12 



(i-«T 

1 + Att + 10t 2 t 2 + 20tH 3 + 10t 4 t 4 + At 5 t 5 + tH 6 

1 + Att + 10t 2 t 2 + 20t 3 t 3 + 35t 4 t 4 + 20tH 5 + 10t 6 t 6 + At 7 t 7 + t 8 t 8 



(1 - tt) 



32 



(2.72) 



The limit N f , N c » 1 and N f = N c + 2. In this limit, (2.67) becomes 



Nt-l 



9 Nf =N c+ 2(t, t, x, y) ~ [1 - [0, . . . , 0, 1; 1, 0, . . . , O^^tt)^ 

+ ([l,0,...,0,l;0,...,0] x . 2/ + [0,...,0;l,0,... ) 0,l] a:;2/ ) (tt) 



Nr 



-[1, 0, . . . , 0; 0, . . . , 0, l] x . y (tt) N f +1 j PE [[1, 0, . . . , 0; 0, . . . , 0, l] x , y tt\ . (2.73) 

Unrefining this Hilbert series, we obtain 

I- N 2 (l-tt) 2 (tt) N J- 1 



9N f ,U(N f -2)(t,t) 



(1 - tt) 



N 2 t 



(2.74) 
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2.4.4 The unrefined Hilbert series for Nf = N c + 3 

Let us focus on the case of Nf = N c + 3. From (2.51), (2.53) and (2.54), we find that 

2 



K Nc (N c ,N c ) 
K Nc (N c ,N c + l) 



(tt) N f + (ttff- 2 (1-tt)' 
+ (tt) N f- l {l-tt) 2 N 2 f , 

-r\t£ft (i-tt)N f 



-r^ttff- 1 (i-tt)' 



N f 
2 



Nf + 1 



ft 



N } 
2 



N f + 1\ 



K Nc {N c , N c + 2) = r 2 (tt) N f (1 - tt) 

K Nc (N c + l,N c + l) = (tt)~ 1+Nf ' 
K Nc (N c + l,N c + 2) = -r\tt) N f (1 - tt) N 
K Nc (N c + 2,N c + 2) = (tt) Nf 



N< 



Nf + 1 



)tt 

2 



tt 



N f 



tt+ (1 -tt) 2 N 2 f 



f > 



and for N c + 2>i>j> N c , the entry K]y c (i,j) is simply K]y c (j,i) with t and t 
interchanged. Other entries of K^ c vanish. 
Thus, it follows that 



det(l - KjvJ = 1- (tt) 



N f -2 



(l-tt)'' 



(tt) 1 



(1 - tt)' 



tt 



2 

Nf + 1 



Nf + 1 



Nf 
2 



-Nj(l - tt) 4 + 2Nj(l - ttftt + 3(tt) 2 



(it)™* . 



tt + tt \2N 2 f (l-tt) + 3tt 

-tt^ 2 +i y 2d -mf 

(2.75) 



Thus, the Hilbert series for the Nf = N c + 3 case is 



9N f ,U{N f -3)(t, t) 



where det(l — fGv f _ 3 ) is given by (2.75). 



det(l - K Nf _, 



(1-tt) 



N f 



(2.76) 
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Examples. We list a few examples of unrefined Hilbert series for small Nf. 

1 + 9tt + 9t 2 t 2 + tH 3 

1 + 9ti + 45t 2 t 2 + 65£ 3 t 3 + 45t 4 t 4 + 9t 5 t 5 + t 6 t 6 
9*, m ~ (— Tie , 

1 + 9tt + 45t 2 t 2 + 165t 3 t 3 + 270t 4 t 4 + 270t 5 t 5 + 165t 6 t 6 + 45t 7 t 7 + 9t 8 t s + t 9 t 9 

96 > m ~ (T^f • 

(2.77) 

The limit Nf, N c » 1 and Nf = N c + 3. To obtain the asymptotic formula, we 
approximate det(l — Kn c ) by neglecting terms of order t 2N f and smaller in comparison 
with those of order t Nf . Expanding the square bracket in the first line of (2.75) and 
neglecting terms of order ( and smaller in comparison with those of order ( ^) , we 
obtain 

1- ( N 2 f) 2 (l-tt) 4 (tt) N f- 2 

a - tty 

2.5 An exact refined Hilbert series for any Nf and N c 

From (2.63) and (2.73), we see that one can express the Hilbert series in terms of a 
sum of representations of SU(Nf) x SU(Nf). We conjecture that the Hilbert series of 
U(N C ) SQCD with Nf flavours can be written as 



9N f ,u(Nf-3){t,t) ~ — ^ . (2. 78) 



X 



g Nf ,u(N c )(t,t,x,y)= (tt) J2a = iana 

ni,ri2,...,njv c =0 

[m, n 2 , ■ ■ ■ , n Nc , 0, . . . , ; 0, . . . , 0, n Nc , . . . , n 2 , n^.y . (2.79) 

A consistency check. We perform a non-trivial check for this formula by setting 
x — y — 1 and obtain an unrefined Hilbert series which can be compared with known 
results. Recall the Weyl dimension formula for the SU(n) irreducible representation: 

dun[ai,...,a n _i] = 1 1 J — . (2.80) 

l<i<j<n ^ % 

For example, let us take Nf = 5 and iV c = 3. Then, 

dim [n u n 2} n 3 ,0;0 } n 3} n 2 ,n 1 ] = (4! 3! 2! l!)- x x 

(m + l)(ni + n 2 + 2)(ni + n 2 + n 3 + 3)(n 1 + n 2 + n 3 + A)x 
(n 2 + l)(n 2 + n 3 + 2)(n 2 + n 3 + + 3)x 
(n 3 + l)(n 3 + + 2)x 

(o + i)l 2 . 
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From (2.79), we find that the unrefined Hilbert series for Nf = 5 and N c = 3 is 
, ~ N 1 + Ml + 10t 2 P + 20t 3 t 3 + 10t 4 t 4 + 4t 5 t 5 + tH & 

95,U(3)(t,t) = -jj . (2.81) 

(1 - tt) 

Observe that this coincides with (2.72). 

The generators and relations. The information about the generators of the moduli 
space and their relations can be extracted from the Hilbert series using the plethystic 
logarithm [3, 4]. To remind the reader, we define the plethystic logarithm of a multi- 
variable function g(t±, ...,t n ) to be 

PL[g(t u . . . , Q) := ^) l 0gg (tl ...,t k n ). (2.82) 

k=l 

where ji{k) is the Mobius function. The significance of the series expansion of the 
plethystic logarithm is stated in [3, 4]: the first terms with plus sign give the generators 
while the first terms with the minus sign give the relations between these generators. 
Now let us compute the plethystic logarithms of various Hilbert series we have computed 
and interpret the results. 

For Nf < N c , we use the Hilbert series is given in the first line of (2.59) and so we 
have 

PL\g Nf < Nc (t,t,x,y)] = [l,0,...,0;0,...,0,l] x;y tt . (2.83) 

This indicates that the moduli space is freely generated by the mesons in the represen- 
tation [1, 0, . . . , 0; 0, . . . , 0, 1] of SU(Nf) x SU(Nf). This agrees with the discussion in 
§2.1. 

For Nf — N c + 1, we use the Hilbert series is given in the first line of (2.63), and 
so we have 

PL[g Nc+l>um {tJ, x, y)\ = [1, 0, . . . , 0; 0, . . . , 0, l] x , y tb - (tt) N < +1 . (2.84) 

This indicates that the moduli space is a complete intersection. The generators are the 
mesons in the representation [1, 0, . . . , 0; 0, . . . , 0, 1] of SU(Nf) x SU(Nf), and there is 
one relation in the trivial representation of SU(Nf) x SU(Nf). This also agrees with 
the discussion in §2.1. 

For Nf > N c + 2, we use the Hilbert series is given in (2.79), and so we have 

PL[g Nf > Nc+2 (t, % x, y)] = [1, 0, . . . , 0; 0, . . . , 0, l) x - y tt 

-[0, . . . , 0, 0, . . . , 0; 0, . . . , 0, l Nc+1 . R , 0, . . . , 0](tt)^ +1 + . . . . (2.85) 
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This indicates that the moduli space is not a complete intersection. The generators are 
the mesons in the representation [1, 0, . . . , 0; 0, . . . , 0, 1] of SU (Nf) x SU(Nf), and there 
are relations in the representation [0, . . . , 0, ljv c +i;i> 0, . . . , 0; 0, . . . , 0, 1jv c +i ; .r, 0, . . . , 0] of 
SU(Nf) x SU(Nf). This also agrees with the discussion in §2.1. 

2.6 Various asymptotics 

In this subsection, we examine asymptotics of Hilbert series for Nf > N c when Nf and 
N c is large. (Note that for Nf < N c , one can simply use the exact formula (2.60).) 

2.6.1 Asymptotics for Nf,N c >> 1 with a fixed difference Nf — N c > 

Let us focus on the limit of large Nf and N c where the difference Nf — N c is kept fixed 
and being positive. For convenience, we define 



We also assume that A > and that A is of order 1. 

In this limit, elements of Kn c are much smaller than 1. Therefore, 

det(l - K Nc ) = exp (Tr log(l - K Nc )) ~ exp (- Tr K Nc ) ~ 1 - Tr K Nc . (2.87) 

Now let us consider TrK^ c . We claim that the leading contribution to TtKn c 
comes from [/jv c +iVL(jv" c +i) in the matrix element Kn c (N c , N c ). Below, we show 
that for < m < A - 1, 



A:— Nf - (N c + 1) . 



(2.86) 



U N c +l+m+\V-(N c +l+m+l) ~ N c U N c +l+mV-(N c +l+m) ■ 



(2.88) 



Note that from (2.51), for m > A + 1, we have UN c +i+ m V-(N c +i+m.) = 0. 
Proof. Suppose that < m < A — 1. Using (2.53) and (2.54). we obtain 



U Nc+1+m +iV HNc+l+m+1) = f [ 2 Fi(m + 1 - A, JV>; N c + 2 + m + 1; tt)] 




x 



Now consider the binomial coefficient 



A — m — 1 




A - m f Nf 



) 



N c + 2 + m VA - m 



Hence, we find that 



UN c +l+m+lV-(N c +l+m+l) 



u A 



N c +l+m,V-(N c +l+m) 



(tt) 



A — m 

N c + 2 + m 



2 F 1 (m + 1 - A, N f ; N c + 2 + m + 1; tt) 
2 F l {m - A, Nf, N c + 2 + m; tt)) 



The ratio between the two hypergeometric functions is of order 1, i.e. as N c , Nf — > oo, 
2 F 1 (m + 1 - A, N f ; N c + 2 + m + 1; tt) 



and so 



2 Fi(m - A, JV /; N c + 2 + m- tt)) 

^AT c +l+m+l^-(Af c +l+m+l) ~ \ A — 171 



Oil) 



(tt) 



+m) 

Since A is of order 1, for a large iV c we have 

UN c +l+m+lV-(N c +l+m+l) 
UN c +l+mV-(N c +l+m) 



N c + 2 + m 



1 



(2.89) 



(2.90) 



(2.91) 



as claimed. Thus, in Tr -ftOv c , we may neglect other diagonal elements of Kn c in compar- 
ison with Kn c (N c , N c ), and the leading contribution in Kn c (N c , N c ) is Un c +iV-^n c +i)- 

□ 

It follows from the above discussion that 

TrK Nc ~K Nc (N c + l,N e +l) 

~ U Nc+ xV-(N c +l) 



( N A" [ 2 F 1 (-A,N f] N c + 2;tt)] 2 (tt) 
l-tt) 2A {tt) N « +1 , 



N c +1 



N f 
A 



(2.92) 



where we have have used the following approximation to establish the fourth equality: 
2 Fx(-A, N f ; N f - A + 1; tt) ~ 2 Fx(-A, JV /; iV /; tt) = (1 - tt) A . (2.93) 
Thus, it follows from the CGBO formula that 



9N f ,U(N c )(t, t) 



i-ra (i-«) 2A («) 



AT C +1 



(1-tt) 



(2.94) 



for Nf,N c » 1 and the difference A = iV / -iV c -l>0 kept fixed. 

Note that for A = 0, we recover the exact formula (2.64) for Nf = N c +1. Moreover, 
for A = 1 and A = 2, we recover the asymptotic formulae (2.74) and (2.78), both of 
which are derived from the exact Hilbert series. 
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2.6.2 Asymptotics for N f ,N c >> 1 with a fixed ratio N f /N c > 1 

In this subsection, we focus on the limit Nf, N c » 1 with a finite ratio 

r:=J->l. (2.95) 

Derivation of the asymptotic formula. Let us consider (2.90). Since A = (r — 
1)N C — 1, we see that for < m < A — 1, we have 

1 A — m , s 

< — < r - 1 . (2.96) 



rN c ~ N c + 2 + m 

Hence, we see that 



tt < U Nc +l+m+lV-(N c +l+m+l) < ~, _ x y ^ ^ 



r 2jy2 ~ f/ A r c+1+m V r _(AT c+ i +m ) 

We would like the upper bound to be sufficiently small, so that we have control over 
UN c +i+ m V-(N c +i+m) for m > 0. We are interested in the following limiting cases: 

• Case 1: Consider r — 1 = o(l) as iV c — > oo. In other words, we take r to be close 
to 1 and tt can take any value between and 1. 

• Case 2: Consider tt = o(l) as iV c — > oo. In other words, we take tt to be small 
and r can take any finite positive value greater than 1. 

In both cases, we see that UN c +i+ m +iV-(N c +i+m+i) can be neglected in comparison with 
U Nc+ i +m V^( Nc+1+m y Thus, the approximation 

det(l - K Nc ) ~ 1 - Tr K Nc ~ 1 - U Nc+1 V_ {Nc+1) (2.98) 

is valid in both limiting cases. We therefore consider 



U Ne+1 V- (Ne+1) = ( N *l^) [2F 1 (N c + l-N f ,N f] N c + 2;tt)Y(tt) 

= (^\)^ 2 Fi(-(r-l)iV c + l,rAr c ;Ar c + 2;tt)] 2 (tt)^ +1 . (2.99) 

Unfortunately, we are not aware of a good asymptotic formula for 
2 Fi(-(r - l)N e + 1, rN c ; N c + 2; tt) 



\N C +1 
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in the limit of large N c , finite r > 1 and < tt < 1. Therefore, we leave the expression 
for the hypergeometric function as it is. 2 
Therefore, we have 

det(l - ~ 1 - (^ + 1 ) a [^r,tt)] a (tt) %+1 , (2-102) 

where ^{N^ r, tt) is a shorthand notation for the hypergeometric function: 

T(N c ,r,t7) :^ 2 F 1 (-(r-l)N c + l,rN c ;N c + 2;t7) . (2.103) 

The asymptotic formula. Thus, using the GCBO formula, we obtain the required 
asymptotic formula 



1 - GSJ 2 \AN^r,tt)] 2 {tt)»«» 

[l-tty 



9rN c ,U(N c) (t,t) VJVC+1; i • (2-104) 



We emphasise that this asymptotic formula is valid for both limiting cases described 
above. We provide non-trivial consistency checks for this formula in the next subsection. 



Special case r = 1: Observe from (2.104) that as r — > 1, we have 

1 

;i ~tt) 



9N c ,u(N c )(t,t) ~ - ^ n2 . (2.105) 



This indeed coincides with the exact result (2.60) for Nf = N c . 



2 On the other hand, the asymptotic formula for the binomial coefficient can be easily computed 
using the Stirling formula logiV! ~ TV log N — N + | log(27riV) for large N, we obtain 

/ rN \ 

log f r , c 1 = log(riV c )! - log(7V c + i)! _ l og ((r - 1)JV C - 1)! 



N n + 1 



N c [r logr - (r - 1) log(r - 1)] + \ log ) ■ l2 l un ) 



Therefore, we have 



rN r 



N c + lJ y/2^N~ c 



(r - l) JV -( r - 1 )- 



(2.101) 



However, the asymptotic formula (2.104) is more elegant when written in terms of the binomial 
coefficient instead of its asymptotic formula. Hence, we leave the binomial coefficient as it is. 
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2.6.3 Consistency checks 

In this subsection, we provide consistency checks for the exact and the asymptotic 
formulae we have derived so far. 

In Figure 1, we plot the graphs of log(?Ar t , + 2,c/(A r c )(^ t) given by (2.74), (2.94) and 
(2.104) against N c , with A = 1, r = 1 + 2/N c and it = 0.5. Note that, for (2.104), this 
is the first limiting case described in §2.6.2, since r — 1 — > as iV c — > oo. It can be see 
that the asymptotic formulae (2.94) and (2.104) approach the exact result (2.74) when 
N c is large. 

l°g(gjV,. + 2,!/(A1,)) 




Figure 1. The case of Nf = N c + 2: The graphs of the exact Hilbert series (2.74) with its 
asymptotic formulae (2.94) and (2.104). We take A = 1, r = 1 + 2/N c , it = 0.5. Note that 
the graphs for both asymptotic formulae are on top of each other in this figure. 

In Figure 2, we plot the graphs of log#jv c +3,£/(jv c )(£> t) given by (2.76), (2.94) and 
(2.104) against iV c , with A = 2, r = l + 3 /N c and it = 0.5. Note that, for (2.104), this 
is the first limiting case described in §2.6.2, since r — 1 — > as N c — > oo. It can be see 
that the asymptotic formulae (2.94) and (2.104) approach the exact result (2.76) when 
N c is large. 

In Figure 3, we plot the graphs of log ge t u(3) (t, t) given by (2.76) and (2.104) (with 
r = 2) against tt. As we expected from the second limiting case in §2.6.2, the asymptotic 
formula should give a good approximation when tt = 0(N~ l ). Indeed, as one can 
see from the graph, the asymptotic result is in agreement with the exact result for 
it < 1/N C = 1/3. 
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10g(gjV„+3,I/(Ay) 




Figure 2. The case of Nf = N c + 3: The graphs of the exact Hilbert series (2.76) with its 
asymptotic formulae (2.94) and (2.104). We take A = 2, r = 1 + 3/N c , tt = 0.5. 



log(&>.(/(3)) 




Figure 3. The case of N c = 3, Nf = 6: The graphs of the exact Hilbert series from (2.76) and 
its asymptotic formula (2.104) (with r = 2). As we expected from the second limiting case in 
§2.6.2, the asymptotic formula is in agreement with the exact result for tt < 1/N C = 1/3. 

3 SU(N C ) SQCD with N f flavours 

Consider M = 1 supersymmetric SU(N C ) gauge theory in four dimensions with Nf 
quark Q\ and Nf anti-quarks Q", where a — 1, . . . , N c and i — 1, . . . , Nf. The super- 
potential of this theory is zero: W = 0. The information about the gauge and global 
symmetries as well as how the matters transform under such symmetries is collected in 
Table 2 . 
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Gauge symmetry 
SU(N C ) 


Global symmetry 

SUiNfr SU(N f ) 2 U(1) B U(1) r U(1)q U(1)q 


Ql 

Qi 


[0,...,0,1] 
[1,0,. ..,0] 


[1,0,. ..,0] [0,...,0] 1 1 
[o,...,o] [0,...,0,1] -1 -1 



U(AMi 




U(N f ) 2 



Table 2. The gauge and global symmetries of SQCD and the quantum numbers of the chiral 
supermultiplets. The quarks are Q\ while the antiquaries are Q®. We also draw it as a quiver 
theory. The circular node represents the U(N C ) gauge symmetry while the two square nodes 
represent global U{Nf)\ and U(Nf)2 symmetries. Each square node gives rise to a baryonic 
U(l) global symmetry, one of which is redundant. We thus have U(1)qq that combine into 
the non-anomalous U(1)b (sum) and anomalous U(1)a (difference). 



SQCD with the SU(N C ) gauge group and Nf flavours is studied and discussed in 
[9, 25-31]. Various geometrical aspects of the moduli space are examined in [9], where 
the Hilbert series are also computed and studied in details. In the following subsection, 
we present a method for computing the Hilbert series and discuss how to rewrite them 
in terms of Toeplitz determinants. We compare the method used in this paper with 
the one used in [9] in §3.1.1. 

3.1 The computations of Hilbert series 

The Hilbert series for SU(N C ) SQCD with Nf flavours can be computed in a similar 
way to that for U (N c ) SQCD with Nf flavours. The first step is to construct the Hilbert 
series of the space of symmetric functions of quarks Q and antiquarks Q. This is given 
by 

N c N c 

t[l,0,...,0] x J^z- 1 + t[0,0,...,l] y J^z a , (3.1) 

a=l a=l 

where for the SU(N C ) gauge group, one needs to impose the condition that 

z 1 z 2 ...z Nc = l. (3.2) 

The second step is to integrate over the gauge group SU(N C ) in order to obtain the 
Hilbert series which counts gauge invariant operators. For this step, we need an ap- 
propriate Haar measure of SU(N C ) for this problem. Let us now discuss such a Haar 
measure. 



PE 
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The Haar measure. The Haar measure of SU (N c ) can be obtained from the Haar 
measure of U(N C ) by restricting Z\ • ■ ■ zn c = 1. This can be written as 



d[lSU(N c ) 



dzi 



N c \(2m) N ^ J lzil=1 Zl /[ , =1 z Ni 



dz Nc 



A Nc (z)\ 2 S( Zl ---z Nc -l) .(3.3) 



In Appendix B, we show that the delta function in this multi-contour integral can be 
replaced by an infinite sum as follows: 3 



8{zx •■•z Nc -l) 



— V z k 
2iri 1 



k „k „fc 



Z . . . Z 



(3.4) 



k=— oo 



Thus, the Haar measure of SU(N C ) can be written as 



/ 



dfi 



SU(N C ) 



N c \{2m 



k=— oo 



dzi 



|*i[=l Z l •/|*W (S [=1 Z W 



A Nc (z)\ 2 4...z k Nc 



oo „ 

/ dVu(N c )4 ■ ■ ■ 



'N c 



(3.5) 



k=—oo 



This formula tells us that the Haar measure of the group SU (N c ) can be written in 
terms of an infinite sum and the U(N C ) Haar measure. 

The Hilbert series. Following from the above discussion, we see that the Hilbert 
series for SU (N c ) SQCD with Nf flavours is given by 



9N f ,su(N c )(t,t,x,y) = / d/j, S u(N c ) PE 



N c 



a=l 



a=l 



k=—oo 



/ dfi u{Nc) z k . . . z k Nc PE 



N c 



N c 



a=l 



0=1 



(3.6) 



Note that the expression in the square bracket is a Toeplitz determinant Dn c (<&) with 
the symbol $ given by 

$(M, x, y, z- k) = z k PE [t[l, 0, . . . , O]^" 1 + t[0, 0, . . . , l] y z] 

= z k <P{t,t,x,y,z) , (3.7) 



3 We are very grateful to Harold Widom for pointing out this substitution to us and providing us 
with a note on analysis of Toeplitz determinant with a non-zero winding symbol. 
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where 4>(t, t, x, y, z) is the symbol for the Toeplitz determinant for U(N C ) SQCD, namely 

0(t,t,x,y,z) = PE[t[l,0,...,0} x z- 1 +t[0,0,...,l} y z] . (3.8) 

Note that, for k ^ 0, $ has a non-zero winding number around the origin, i.e. 

y 1 f dz d z $(t,t,x,y,z;k) 
2m J z $(t,t,x,y,z;k) 

Note however that the method we have used in §2.2.1 applies only to symbols with zero 
winding number. Therefore, such a method needs to be generalised in order to compute 
the Toeplitz determinant with the symbol $. In the next subsection, we discuss related 
theorems for the cases of non-zero winding numbers. For now, we rewrite the Hilbert 
series of SU(N C ) SQCD with Nf flavours as 

oo 

9N f ,su(N c )(t,t,x,y) = D Nc ((f)) + Y^ [ D N c (z k <t>(t,t,x,y,z)) + D Nc (z~ k (j)(t,t,x } y, z))] 

k=i 

oo 

= 9N f ,u(N c )(tJ,x,y) + [ D N c (z k <l>(t,t,x,y } z)) + D Nc (z~ k <t>(t,t,x,y,z))] .(3.10) 

fc=i 

This tells us that the Hilbert series of SU(N C ) SQCD with Nf flavours is simply the 
Hilbert series of U(N C ) SQCD with Nf flavours with correction terms coming from 
non-zero winding parts of the symbol. Subsequently, we refer to the first, second and 
third terms in (3.10) respectively as the zero winding part, the positive winding 
part and the negative winding part. 

The unrefined Hilbert series. Setting x\, . . . ,xn c and yi, . . . , y^ c to unity, we have 
the unrefined Hilbert series: 



9N f ,SU(N c )(t,t) = 



k=— oo 



Za 



(3.11) 



Note that the expression in the square bracket is a Toeplitz determinant Dn c (§) with 
the symbol $ given by 



z k 



$(t,t,z;k) = z k <p(t,t } z) = ^— — — . (3.12) 

v ' ) ^ ' ) (l-tz) N f(l - ^) N f 
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3.1.1 Comparison with the method used in the 'Apergu paper' 

We emphasise that the method of computation we use in this paper is significantly 
different from the one in the 'Apergu paper' [9]. In this paper, we use coordinates 
z Nc on the maximal torus of U(N C ) and then impose the condition Yl a =i Za = 1 
using a delta function. In this way, one can recast the Molien-Weyl formula into a 
Toeplitz determinant, which is a key tool for computations. There are a number of 
techniques from the random matrix theory that can be used to evaluate such Toeplitz 
determinants both exactly and asymptotically for a large class of values of N c and Nf. 
We present these techniques and apply them to our computations in the following 
subsections. 

In [9], on the other hand, the integrand and the SU(N C ) Haar measure are written 
in terms of N c — 1 coordinates on the maximal torus of SU(N C ) without a delta function. 
In this way, one can compute Hilbert series case by case for a given (Nf,N c ). These 
computations become very cumbersome when N c is large due to a large number of 
contour integrals. However, based on case by case results, one can make conjectures 
about the general results, and indeed in [9] a number of these are stated as observations. 
In the following subsections, a number of these observations can be proven (or at least 
can be checked in a non-trivial way) using Toeplitz determinants. Moreover, we show 
that this technique can be used to compute a number of new exact Hilbert series and 
asymptotic formulae for large Nf and N c - many of which are too difficult or impossible 
to be derived using the method in [9]. 

3.2 The Bottcher-Widom and the Fisher-Hartwig theorems 

In this section, we state two theorems which are very useful in computing the Toeplitz 
determinant with a symbol which has non-zero winding number around the origin. In 
what follows, we still use the same definitions of the Toeplitz matrix, the Hankel matrix, 
U, V, G(4>) and E[<f>) as in §§2.2.1 and 2.2.2. 

Theorem 3.3. (The Bottcher-Widom theorem. [40]) Let : T — > C — {0} be a 

continuous function such that <p(z) ^ for all z e T and <p has a zero winding number 
around the origin, i.e. 




(3.13) 



• If k > and T n+ k(4>) is invertible, then both of the operators 



1 - H(b)H(Z)Q n+k and 1 - H(b)Q n H(c)Q k 



(3.14) 



30 



are invertible, and 

D n (z- k <p) = {-l) nk D n+k {<j>)F n>k {4>) , (3.15) 

where 

F„, fc (0) = detP fc (l - H{U)Q n H{V)Q k )- l T{z- n U)P k . (3.16) 

• On the other hand, assuming k > 0, we have 

D n (z k <j>) = (-l) nk D n+k {<f>)G(<f>)- k G(U) k F n> M , (3.17) 

where 

F n ,S) = detP k T(z n V)(l - Q k H{U)Q n H{V)Y 1 P k . (3.18) 

The following theorem gives the asymptotic formula for D n {z~ k (f) when n is large. 

Theorem 3.4. (The Fisher-Hartwig(-Bottcher-Silbermann) theorem. [46— 
48]) Assume <p to be as in the Bottcher-Widom Theorem. 

• If k > 0, then 

F n , fc (0) ~ detT k (z- n U) . (3.19) 

It then follows from the Bottcher-Widom theorem and the strong Szego limit 
theorem that 

D n (z- k <p) ~ (-l) nk G(<P) n+k E(<p) det T k {z~ n U) . (3.20) 

• On the other hand, assuming k > 0, we have 

F ntk ($) ~ detT k (z n V) , (3.21) 

It then follows from the Bottcher-Widom theorem and the strong Szego limit 
theorem that 

D n (z k <P) ~ {-l) nk G{<j)) n E{<f))G{U) k detT k {z n V) . (3.22) 
Subsequently, we apply these theorems to compute Hilbert series of SQCD. 
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3.2.1 Applications to the Hilbert series of SQCD 

Recall the discussion around (3.7). We take to be the symbol that we have used for 
U(N C ) SQCD, i.e. 

(j)(t, t, x, y, z) = PE [t[l, 0, . . . , O]^" 1 + t[0, 0, . . . , l) y z] . 

Then, the symbol $ for the Toeplitz determinant of SU(N C ) SQCD is 

$(t, t, x, y, z; k) = z k PE [t[l, 0, . . . , O]^" 1 + t[Q, 0, . . . , l] y z] 
= z k (j)(t,t,x,y,z) . 

Hence, we can apply the Bottcher-Widom and Fisher-Hartwig theorems to compute 
the N c x N c Toeplitz determinant with the symbol t, x, y, z\ k). Thus, we take n in 
these theorems to be the number of colours N c . 

Results for any Nf and any A^ c . With the symbol $(t,t,x,y, z;k), we find that 
for any N c > 1, Nf > 1 and any k G Z, 

Q Nc H(V)Q k = , 
Q k H(U)Q Nc = , 

G(U) = 1 . (3.23) 

Let us now look at various cases of Nf and N c . 

3.2.2 The case of JV> < N c - 1 

For Nf < N c — 1, we find that, for all k > 0, 

P k T(z- Nc U)P k = , (3.24) 

and so it follows from the Bottcher-Widom theorem that 

F Nc , k {4>) = . (3.25) 

In other words, the Toeplitz determinants vanish for all negative winding number of 
the symbol $, i.e. for all k > 0, 

D Nc ($(t,t,x,y,z;-k)) = D Nc ( Z - k (j>(t,t,x,y,z)) = . (3.26) 

Similarly, for all k > 0, we have 

F Nc ,S) = , (3.27) 
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and hence the Toeplitz determinants vanish for all positive winding number of the 
symbol i.e. for all k > 0, 

D Nc ($(t, t, x, y, z-k))= D Nc (z k <P(t, t,x,y,z))=0. (3.28) 

Therefore, it follows from (3.10) that the Hilbert series for SU(N C ) SQCD with 
Nf < N c — 1 flavours is exactly equal to the Hilbert series for U(N C ) SQCD with 
N f < N c - 1 flavours: 



9N f < Nc -i(t, t, x, y) = PE [[1, 0, . . . , 0; 0, . . . , 0, l] x , y tt\ 

22 [ni,n2,...,n Nf -i;n Nf ^i,...,n2,n 1 ](tt)^= 1 Jrtj '. (3.29) 



1 - (tt) N f ^ n 

The unrefined Hilbert series is then 

g Nf <N c -i(t,t) = - — -—s . (3.30) 
(1 - tt) N f 

Note that these Hilbert series are in agreement with (3.6) and (5.3) of [9]. 

The moduli space for Nf < N c — 1. As can be seen from the Hilbert series, the 
moduli space of SU(N C ) SQCD with Nf < N c — 1 is freely generated by the meson 
transforming in the bi-fundamental representation [1, 0, . . . , 0; 0, . . . , 0, 1] of SU(Nf) x 
SU(Nf). 

Let us make a brief comment on quantum corrections. A non-perturbative Aflieck- 
Dine-Seiberg (ADS) superpotential [25, 28, 29, 32] is dynamically generated. This 
completely lifts the vacuum degeneracy and hence there is no supersymmetric vacuum. 
While the Hilbert series of the classical moduli space does not have a physical meaning 
in the full quantum theory, it nevertheless contains information about the structure of 
gauge invariant operators for Nf < N c . 

3.2.3 The case of N f = N c 

It follows from the Bottcher-Widom theorem that 

F Nc , k (<P) = (~t) kNc , F Nc , k {4>) = (~t) kNc ■ (3-31) 
Then, from (3.15), for k > 0, 

D Nc {z- k <P) = (-l) kN «D Nc+k (4>)(-t) m < = Av c+fc (# % . (3.32) 

Recall that for U(N C + k) SQCD with N f = N c flavours, 

D N c +k(4>) = 9N c ,u(N c +k)(t,t,x,y) =PE [[1,0,...,0;0,...,0, l] x . y tt\ . (3.33) 
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Thus, we have 

oo 

V D No (z- k <f>) = PE [[1, 0, . . . , 0; 0, . . . , 0, l] x , y tt\ . (3.34) 

, — : 1 — t c 



k=l 



Similarly, using (3.17), we find that 

E D *° = Y~t^ PE °' • • • ' 0; °' • • • ' °' 1]x ' yt ^ ■ (3 - 35) 

fc=l 

The Hilbert series for SU(N C ) SQCD with Nf = N c flavours is then given by (3.10): 

/ t Nc t Nc \ 

g Nf ,su(N c) (t,t,x,y) = \ 1 + + JN* ) PE[[l,0,...,0;0,...,0,l] xiW tt] 



1 _ (u) Nc 

PE[[l,0,...,0;0,...,0,l} x . y tt\ 



(l-t N *)(l -t N *) 

= (1 - {tt) N °) PE [[1, 0, . . . , 0; 0, . . . , 0, l] x . y tt + t N ° + P c ] . (3.36) 

This can be written in terms of a sum of irreducible representations of SU (Nf) x SU (Nf) 
as 

9N c ,SU(N c ){t,t,X,y) = 2^ * ' c t^=i J 

ni,ri2,...,njv c -i,^,m.>0 

[m, n 2 , • • • , njv c -i ; njv c -i, ■ ■■ ,n 2 , ni] x;y . (3.37) 

This is in agreement with the general formula (5.2) of [9]. 

Setting x's and y's to unity, we obtain the unrefined Hilbert series: 

1 _ (U)N C 

9n su(N)(t,t,l,l) = ^—o — - — • (3-38) 

This is in agreement with (3.15) of [9]. 

The moduli space for Nf = N c . Indeed, this Hilbert series (3.36) tells us that the 
moduli space is a complete intersection (see, e.g., [9]). The generators are mesons 

M) = QIQ" , (3.39) 

transforming in the bi-fundamental representation [1, 0, . . . , 0; 0, . . . , 0, 1] of SU (Nf) x 
SU(Nf), and baryons and anti-baryons 

B n.., Nc = e au..a NoQ n . . . Q m j = ^^Q* . . . QjJ- , ( 3 .40) 
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each of which transforms as a singlet under SU(Nf) x SU(Nf). There is one relation 
(at order (tt) Nc ) between these generators, namely 

detM- (*B)(*B) = , (3.41) 

where *B = -^-,ej 1 ...j JVc -B n " : ' JVc and similarly for *B. Observe that the baryons, anti- 
baryons and their relation with mesons come from correction terms to the Hilbert series 
of U (N c ) gauge theories with Nj < N c flavours. 

Now let us make a brief comment on the quantum moduli space. It is still generated 
by M, B and B. However, the classical relations (3.41) is modified by a one instanton 
effect [25-31], and the quantum moduli space is described by the relation 

det M - {*B) (*B) = A 2iVc , (3.42) 

where A is the scale of the theory. Although details of the relation are modified, the 
representations in which the generators and their relation transform under the global 
symmetry are unaffected. Thus, in spite of different geometrical properties between 
the classical and quantum moduli spaces, the Hilbert series is not corrected quantum 
mechanically. 

3.2.4 The case of N f = N c + 1 

It follows from the Bottcher-Widom theorem that 

k 

F Nc ,M = (~t) kNc [°> • • • ' °> W> °> • • • ' 5 k-m,0,..., 0] x , y (-tt) m , 

m=0 
k 

F Nc ,S) = {-t) kN < £)[0, 0, . . . , k - m ; 0, . . . , 0, l m , R , 0, . . . , 0} x , y {-tt) m . (3.43) 

m=0 

For k > 0, the Toeplitz determinant D Nc+k ((p) is given by 

D Nc+k (<i>)=g No+imNo+k) (tj,x,y) = PE[[l,0^ , (3.44) 

and the Toeplitz determinant Dn c (4>) is given by 

DnM) = 9N c+ i,u(N c )(tX x, y) = (1 - (tt)^ +1 ) PE [[1, 0, . . . , 0; 0, . . . , 0, l] x , y tt\ (3.45) 

Thus, from (3.10), we find that the Hilbert series of SU(N C ) SQCD with N f = N c + 1 
flavours can be written as 



9N c +i,su(N c )(t,t,x,y) 



1 _ (ttfc+i + >T J2 [°. ■ • • > °> W« 0) ■ • • >0 ; k — m,0, . . . , 0] BlB (-tt) f 



Y;t kNc J2[0,...,0,k-m; 0,...,0,l m . R ,0,...,0] Xj y(-tt) 

k=l m=0 



k=l m=0 

PE[[l,0,...,0;0,...,0,l] x;y tt] . 

(3.46) 
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This can be rewritten in terms of a sum of irreducible representations of SU(Nf) x 
SU(N f ) as 



" V 1 jrij+mNc x 



9n c+ i,su(n c ) (t, t, x,y)= Yl jnj+m ° ^ jT 

n 1 ,n 2 ,...,n Nc -i,l,m>0 

[m, n 2 ,..., n Nc -x, £ ; m, n Nc -i, ...,n 2 , ni\ X]y . (3.47) 
This is in agreement with the general formula (5.26) of [9]. 
The unrefined Hilbert series. Setting x's and y's to unity, we find that 



' 'N f \ fN f + k-m-l 



k — in 



-tty 



m=0 ^ ' ^ 

= (-*)»• ( Nf + * ~ M M-k, -N,; 1-k-N,; it) , 

^)=H)-E©r /+ **-" m m_i ) ( -«T 

m=0 V / \ / 

= (-t) kN ^ Nf+ ^~ i yF l (-k,~N f] l-k-N f ,tt) . (3.4* 
Then, we find that the negative winding part is 

E ^„(^V) = £^ f ^ + * ~ * W"*' - JV>; 1-A:-JV>; tt) 

k=i I 1 " tt ) ' fc=i V fc y 

(l-tti^Cl-^^IS W + n + V 



(3.49) 



whereas the positive winding part is 



fe=i 

7V C W> 



m=0 n=0 

and from (2.64), we find that the zero winding part is 



c 1 / N \ 

m — n r,— n \777. "T 71 T" / 



1 - (tt) N f 
{l-ttff 



DnM = - %p ■ (3.51) 
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Thus, from (3.10), we obtain the unrefined Hilbert series for SU(N C ) with Nf = N c + 1 
flavours as 



( T\ 1 

9N c +l,SU(N c ){t,t) 



1 - (tt) 



N c +1 



m=0 n=0 



(1 _ tt )(JV c +l) 2 

, N c + 1 \ f t( m+n+1 ') Nc f~(m+n+l)N c 
^-«)2^ Z^(- 1 ) m+ "l m+ C n + lJ ( (1_^ C )JV C +1 + (1 _£v c) JV c+ l ) ^ 



(3.52) 

Examples. For simplicity, we set t = t. From (3.52), we have 
/ x 14" 6t 2 + 6t 4 + t 6 

03,51/(2)1* J - — g , 

(1 -t 2 ) 

1 + 4t 2 + 4t 3 + 10t 4 + 8t 5 + 14t 6 + 8f + 10t 8 + 4t 9 + 4t 10 + t 12 

94,SU(3){t) - (1 -t 2 ) 12 (l -t3)4 ' 

1 + 5t 2 + 20t 4 + 50t 6 + 85t 8 + 100t 10 + 85t 12 + 50t 14 + 20t 16 + 5t 18 + t 20 



95,SU(4)(t) 



i - 1 2 ) 25 {1 + t 2 y 



Note that the first two results are in agreement with those in [9]. These examples 
provide a consistency check of the formula (3.52). We emphasise here that with the 
knowledge of Toeplitz matrices and their determinants, one can obtain exact results 
without encountering difficulties from computing a large number of contour integrals. 

The moduli space for Nf > N c . For completeness of the paper, let us briefly 
summarise the information about the moduli space for Nf > N c (see also [25-31]). 
This information is contained in the Hilbert series and can be extracted using the 
plethystic logarithm. Since this has already been shown in [9], we simply state the 
results here. The generators are mesons 

M) = QiQ« , (3.53) 

transforming in the bi-fundamental representation [1, 0, . . . , 0; 0, . . . , 0, 1] of SU {Nf) x 
SU(Nf), and baryons and anti-baryons 

B n...i Nc = e ai...a NcQ n ? ^ ^ = ^ ^ . (3 54) 

transforming respectively in the representation 

[0,. .. ,0,1^,0,. ..,0;0,...,0] and [0, . . . , 0; 0, . . . , 0, l NeiR , 0, . . . , 0] (3.55) 
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of SU(Nf) x SU(Nf). They are subject to the relations: 

{*B)B = *(M Nc ) , M-*B = 0, M-*B = 0. (3.56) 

where (*B) iffa+1 _ iN = -^e^...^. -B Jl '" lJVc and a '•' denotes a contraction of an upper with 
a lower flavour index. These relations transform respectively in the representations 

[0,...,0,l^ c;L ,0,...,0;0,...,0,l^,0,...,0] , 

[0, 0, l ( 7v c+ i);L, 0, 0; 0, 0, 1] , [1, 0, 0; 0, 0, l {Nc+l);R , 0, 0] . (3.57) 

Note that, in this case, the quantum moduli space coincides with the classical moduli 
space [25, 30]. Thus, the Hilbert series of the classical moduli space for Nf > N c is also 
valid for the quantum moduli space of the theory. 

3.3 An exact refined Hilbert series for any Nf and N c 

From (3.29), (3.37) and (3.47), one can see that the results we have obtained so far are 
in agreement with the general formula (5.26) of [9] which gives the Hilbert series for 
any Nf and any iV c : 

9N f ,su m (t, t) = ^ jni+mc ^ Jnj+mNc 

ni,n2,...,n^r c _i,£,m>0 

[m, n 2 , • • • , n Nc -x, £ Nc;L , 0, . . . , 0; 0, . . . , 0, m Nc . R , njv c -i, • • • , n 2 , ni] . (3.58) 

So far we have discuss a number of exact results for SU(N C ) SQCD with Nf flavours. 
In the next subsection, we study asymptotics of the Hilbert series for large Nf and iV c . 

3.4 Asymptotics for large Nf and N c 

In the limit of large Nf and iV c , it is convenient to obtain the Hilbert series using the 
Fisher-Hartwig Theorem 3.4. 

The negative winding part. Let us first consider the negative winding part in 
(3.10). We would like to evaluate the sum Y^k=i -^ ) A r c( 2;_fc 0) using (3.20). (Recall that 
G{4>) = 1 and G(U) — 1.) We claim that the leading contribution to this sum comes 
from the term with k — 1, namely 

D N £z- x 4>) = (-1)^(0) detT^z-^U) = (-1)^(0)^(^17) . 

This is because for Nf,N c » 1, the leading behaviour of det Tk(z~ Nc U) is t^ Nc ; 
therefore D^ c {z~ k (()) is of order t kNc E((j)) J and hence the other terms in the sum 
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X^fcLi Dn c {z k 4>) can be neglected in comparison with D^ c (z It follows imme- 

diately from the definition (2.21) of the Toeplitz matrix that 

T x {z- N 'U) = ^-l -z~ N 'U = U Nc . (3.59) 

where we have used (2.27) to establish the last equality. Therefore, using (3.20) and 
recalling from (2.32) that 

£(0) = PE[[l,O,...,O;O,...,O,l]^tt] , 

we obtain 

oo 

Y,D Nc {z- k ct>) ~D Nb {z- x <I>) 

k=l 

= (-lf^E^T^z-^U) 

= (-l) N °U Nc PE [[1, 0, . . . , 0; 0, . . . , 0, l\ x , y tt\ , (3.60) 
where we have neglected terms of order 0(t c )E(<j)) and smaller. 

The positive winding part. A similar argument can be applied to compute the 
positive winding part of (3.10). From (2.21) and (2.27), we find that 

detTAz^V) = TAz^V) = — I —z N <V = V- N . (3.61) 

^J\z\ = l Z 

Therefore, from (3.22), we obtain 

oo 

^2 D N c {z k (f>) ~ D Nc {z4>) 

k=l 

= {-l) N °E{<l>)T x (z N *V) 

= (-1) N °V- Ne PE [[1, 0, . . . , 0; 0, . . . , 0, l] x;y tt\ , (3.62) 
where we have neglected terms of order 0(t 2Nc )E((j)) and smaller. 

The zero winding part. The zero winding part Dn c {4>) is discussed in §2.6. Recall 
that there are 3 limiting cases of our interest: 

• The difference A := N f - (N c + 1) > is finite when Nf, N c — >■ oo. 

• The ratio r := Nf/N c = 1 + o(l) as iV c — > oo. (This is case 1 in §2.6.2.) 

• The ratio r > 1 is finite and tt = o(l) as N c — > oo. (This is case 2 in §2.6.2.) 
In these limiting cases, the following approximation is valid: 

D JVc (0)~(l-£/^ + iy_ ( ^ e+1) )PE[[l,O,...,O;O,...,O,l] ! « ; yttl • (3.63) 
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The asymptotic formula. Thus, from (3.10), we find that 

g Nf ,su(N c ){tJ,x,y) ~ (l + {-l) Nc U Nc + {-l) Nc V_ Nc - U Nc+1 V HNc+1) ) x 

PE[[l,0,...,0;0,...,0,l] x;? ,tt] . (3.64) 

3.4.1 Asymptotics for Nf,N c » 1 with a fixed difference Nf — N c > 

Now let us focus on the limit Nf ,N C »1 with a fixed difference 

5 ■= N f - N c > . (3.65) 
From (2.53) and (2.54), we have the following asymptotic formulae: 



Un c = (-t) JVc ( )2Fi(-S,Nf,N f -6 + l;U) 



V- Na ~(^yi-U)'(-tf'. (3.66) 

Substituting (3.66) and (2.92) into (3.64), we find that the asymptotic formula for the 
unrefined Hilbert series is 

t ~\ 1 

9N c +S,SU(N c )(t,t) ~ (1 _ tf)(JVc+(5)2 >< 

1 + ( Nc + 5 )(1 -tt) s (t N " + t N ^ - f^ c + f)V-tt) 2(5+1 W c+1 • (3-67) 

V * *c / V c "I - 1/ 

Special case of = N c . In this special case, the formula (3.67) reduces to 



9n c su(n c ) {t, t, x, y) ^~ jy2 . (3.68) 



1 + t"* + t N ' 
(1 - «) 

It is worthwhile comparing this with (3.38). We see that in the limit of large N c , 



1 - {tt) 



1 + t Nc + t Nc . (3.69) 



(1 -t N *)(l -t N -) 



Substituting this into (3.38), we recover (3.68) as expected. 
We provide other consistency checks for (3.67) in §3.4.3. 
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3.4.2 Asymptotics for N f ,N c » 1 with a fixed ratio N f /N c > 1 

Let r be the the ratio Nf/N c between Nf and N c and assume that r > 1. We apply 
(3.64) to find the asymptotic formula. 
Recall from (2.99) and (2.101) that 



rN, x 2 



i) 



Un M V hnm) ~ [ N '\\ ) [ 3 F!(-(r - l)N c + l,rN c ;N c + 2;tt)] 2 (tt) NM .(3.70) 



where ^(Nc, r, tt) is defined as 

JF(iV c , r, tt) := 2 F x (-(r - 1)JV C + 1, rN c ; N c + 2; tt) . 

Now let us compute (— 1) c [/jv c + (— 1) ^_jv c - From (2.53) and (2.54), we can write 

~ frN \ ~ 
Un c = H) Nc ( N c ) 2 F 1 (-(r - 1)N C , rN c ; N c + 1; tt) , 

V- Nc = {-t) N °[ " ) 2 F 1 (-(r-l)N c ,rN c ;N c +l;tt) . (3.71) 



Writing 



we have 



F(N c ,r,tt) ■.= 2 F 1 (-(r-l)N c ,rN c ;N c +l;tt) , (3.72) 



{-l) N «U Nc + (-l) N <V_ Nc ~ ^ c Y(N c ,r,tt) (t^ + t N ^) . (3.73) 

Thus, from (3.64), we arrive at the asymptotic formula 

1 



9rN c ,SU(N c )(t,t) 



l+( r ^-F(iV c ,r,tt)(t^+^) 



(i - tty 2N c 



N c +l 



(3.74) 



One interesting application of this asymptotic formula is that one can use it to study 
the moduli space in the conformal window, namely for |iV c < Nf < 3N C , where the 
theory has a non-trivial IR fixed point. In the conformal window, the theory possesses 
a dual description, known as the Seiberg duality [27]. It is an interesting problem to 
check this duality using Hilbert series. So far Romelsberger [8] has showed that the 
Hilbert series of SU (2) SQCD with 3 flavours and its magnetic dual match. Now that 
the Hilbert series for the theory in the conformal window is available in the limit of 
large iV c and Nf, the checking should be possible for such an asymptotic limit. Since 
this problem deserves investigation in its own right, we leave this to future work. 
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3.4.3 Consistency checks 

In this subsection, we provide consistency checks for various asymptotic formulae we 
have derived so far. 

In Figure 4, we plot the graphs of loggN c +i,su(N c )(t,t) given by (3.52), (3.67) and 
(3.74) against N c , with 5 = 1, r = 1 + 1/N C , t = 0.3 and t — 0.1. It can be see that 
the asymptotic formulae (3.67) and (3.74) approach the exact result (3.52) when N c is 
large. 

logCftv^+i.sufAy) 




0.4 - 

0.2 - 



Figure 4. The case of Nf = N c + 1: The graphs of log57v c +i,S(7(Af c )(^ t) given by (3.52) with 
its asymptotic formulae (3.67) and (3.74). We take 5 = 1, r = 1 + l/iV c , t = 0.3 and t = 0.1. 
Note that the graphs for both asymptotic formulae are on top of each other in this figure. 

In Figure 5, we plot the graphs of log <7jv c +5,s[/(iVc)(^ t) given by the asymptotic 
formulae (3.67) and (3.74) against iV c , with 5 = 5, r = 1 + 5/N c , t = 0.3 and t = 0.1. 
It can be see that the asymptotic formula (3.67) approaches the asymptotic formula 
(3.74) for large N c . 

Let us now set t = t. In Figure 6, we plot the graphs of log ge t su(3) (t, t) given by 
(3.58) and (3.74) (with r = 2 and N c = 3) against t. As we expected from the second 
limiting case in §2.6.2, the asymptotic formula should give a good approximation when 

1 /2 

t = 0(Nc 1,A ). Indeed, as one can see from the graph, the asymptotic result is in 
agreement with the exact result for t < l/y/N~ c = 1/a/3 ~ 0.58. 
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l0g(giv c +5,SU(JV c )) 




l - 



Figure 5. The case of Nj = N c + 5: The graphs of logg Nc+5 SU ^ Nc ^(t, t) given by the 
asymptotic formulae (3.67) and (3.74). We take 5 = 5, r = 1 + 5/N c , t = 0.3 and t = 0.1. 

log(^6,SU (3)) 




Figure 6. The case of N c = 3,iV/ = 6 and t = t: The graphs of log 56,51/(3) (*> *) horn 
(3.58) and its asymptotic formula (3.74) (with r = 2 and _/V c = 3). As we expected from the 
second limiting case in §2.6.2, the asymptotic formula is in agreement with the exact result 
for t < l/VN~ c = l/>/3 0.58. 
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A Multi-contour integrals and determinants 

In this section, we show that one can rewrite the contour integrals in the Hilbert series 
of U(N C ) SQCD in terms of a determinant. A key tool is Gram's formula (see, e.g., 
Appendix A. 12 of [42]), which states as follows. Let 

n 

b ij = ^2 v ^v* ja , i, j = 1,2, ...,m ; (A.l) 

a=l 

then 

j n n 

— ^ • • • |det(u iaj ) 1 <i !i < m | = det(foj i )i< iiJ < m . (A.2) 

ai=l a m =l 

The U (N c ) Haar measure. Let us first consider the U(N C ) Haar measure 
^u m = Ncl( l i)Nc { ^|A. C (,)' 2 



[«l|=l z l J\z Nc \ = l Z N, 



i I dzi I &Z Nc | b _ 1 ,2 , . . 

det{z a )i< ,6<jvj • (A.3) 



N c \(2m) N - J {zil=1 zi J\z Nc \=i zn, 
Let us now apply (A.2). We take m = N c and make the following substitutions: 



E 

a=l 
n n 

E-E 

ai=l o m =l 



i r dz 



2-Ki ./j z | =1 z 



[2m) N ° J\ zi]=1 z x J\ ZNc \=\ z Nc 



v iaj — > z*- 1 . (A.4) 
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Thus, from (A.l), we find that 

b ab = ^-J ^-\z b -y = ±<f ^"V" 6 = -L / ^z«- b . (A.5) 
2m J lzl=1 z 2m J N=1 z 2m J ]z]=1 z 

Thus, from (A. 2), (A. 3) and (A.5), we can rewrite the Haar measure of U(N C ) as 

d/i u(Nc) = det ( / . (A.6) 

The [/ (N c ) SQCD Hilbert series. Now consider the integral of the type: 

11^^)" (A-7) 

a=l 

Note that this is also the integral we encounter in (2.35) in the context of U(N C ) SQCD. 
We make a similar substitution as in (A. 4): 

y — ► — / —<KM 



n n - „ , „ , N c 

— ► ( A - 8 ) 

Thus, from (A. 2), we find that 

AO = det (-L t ^z a ~ b ^,z)) . (A.9) 



27Ci Az\=l Z / l <a ,b<N r 



B A delta function in the contour integral 

Let f(z) a function which is analytic everywhere on the unit circle. In this section, we 
show that the delta function in the integral 



'1=1 

can be replaced by an infinite sum as 4 



I ^ § (z-l)f(z) (B.l) 

J\z\=l Z 



1 oo 



2m 

k=—oo 



1 We are very grateful to Harold Widom for pointing out this substitution to us. 
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so that we have 



I |^-i)/(z) = 5 L f; / (b.3) 

J\z\=l Z lm J\z\=l Z 



'1-1 = 1 * k= _^J\z\ 

Proof. Consider 

= I ^ 6 (z _ l)f{z) =i r d0 8{d e - l)f(e^) . (B.4) 

J\Z\ = 1 Z J-7T 

Note that 

oo . oo 



2tt 

n=— oo fc=— oo 



where the first equality can be verified by integrating from l-nm — n to 27rm + 7r (for 
any m € Z) and noting that —i = J_ d9 S(e ld — 1), whereas in the last step we have 
used the Poisson summation formula. Substituting this back into (B.4), we obtain 

I ^ 5{z _ l)f{z) = 1_ V r d9 e ik6 f(e id ) . (B.6) 
Writing z = e ld on the right hand side, we arrive at 

This amounts to the replacement 

1 oo 
fc=— oo 
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